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Abstract 



We calculate the lower moments of the deep-inelastic structure 
functions of the pion and the rho meson on the lattice. Of particular 
interest to us are the spin-dependent structure functions of the rho. 
The calculations are done with Wilson fermions and for three values 
of the quark mass, so that we can perform an extrapolation to the 
chiral limit. 

1 Introduction 

The deep-inelastic structure functions of nucleons and mesons are currently 
receiving a lot of attention, both experimentally and theoretically. It has 
become possible to compute these structure functions from first principles 
within the framework of lattice QCD [1, 2]. This allows a quantitative test 
of QCD which goes beyond perturbation theory. The basis of the calculation 
is the operator product expansion which relates the moments of the struc- 
ture functions to forward hadron matrix elements of certain local operators. 
Lattice simulations of these matrix elements, combined with an appropriate 
calculation of the Wilson coefficients, can in principle provide complete in- 
formation of the quark and gluonic structure of the hadronic states. The 
aim of this paper is to compute the structure functions of the pion and the 
rho meson. Following [2] we use Wilson fermions and work in the quenched 
approximation, where internal quark loops are neglected. 

The pion structure function, which so far was extracted from measure- 
ments of the Drell-Yan lepton-pair production cross section [3], is directly 
being measured at HERA at present [4], and we expect that new data will 
become available soon. Using the techniques described in [5, 6], it should in 
principle be possible to measure the structure functions of the rho as well in 
the near future. Some of the spin-dependent structure functions, in particu- 
lar, should be easy to separate from the dominant pion exchange process. 

But even without having any experimental data to compare with, the 
internal structure of the rho meson is an interesting subject to study on the 
lattice. Besides the structure functions already known from the nucleon, 
one finds new structure functions that contain qualitatively new information 
which has no analog in the case of spin-i targets [7]. We hope that this 



2 



investigation will lead to a better understanding of quark binding effects in 
hadrons. 

The structure functions of the rho are also of interest for the interpretation 
of photoproduction and two-photon inclusive cross sections, as the photon 
has a substantial hadronic component which to a good approximation can 
be described by the rho meson. 

For a spin-0 target like the pion, the kinematical framework is simpler 
than in the familiar nucleon case. The details for the case of a polarized spin- 
1 particle have been worked out by Hoodbhoy et al. [7] . The hadronic tensor, 
i.e. the imaginary part of the forward current-hadron scattering amplitude, 

W^(p,q,X,X') = ±-J d 4 xe^ x (p,A'|^(^),r(0)]b,A) (1) 

(with A, A' labeling the polarization) decomposes into eight structure func- 
tions: 

J) D 1 I 

+ \b A {s^ - t^) + i—t^ X(T q x s a 

+ i^x«q\p-qs° - s ■ qp°) , (2) 

where v — p ■ q, and r^ u , s^ u , t^ u , u^ v are kinematical tensors [7] constructed 
from the momentum transfer q and the polarization vector e. The latter 
satisfies e • p = 0, e 2 = — m 2 , and 

S<T = ^ a " T ^ , (3) 

with m being the hadron mass. Here ^ vXa is the completely antisymmetric 
tensor with e 0123 = — 1. 

The structure functions F±, F 2 , g±, and g 2 play the same role as for a spin-| 
target. In the parton model, the structure of the hadron can be described by 
the probability q™{x) (q™(x)) of finding a quark with momentum fraction x 
and spin up (down) along the direction of motion when the hadron is moving 
with infinite momentum and has spin projection m = 0, ±1. Symmetry 
implies 

qi(x) = ql\x), q\(x) = q^(x) , 

<#(*) = , (4) 
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so there remain only three independent parton distribution functions ql(x), 

q\(x), and q®(x). In leading order, the single-flavor structure function Fi 9 \x) 
is one half of the probability to find a quark q with momentum fraction x, 
and F 2 q \x) obeys the Callan-Gross relation: 

Fi q \x) = ±(q*(x) + ql(z) + <${x 
F 2 [q \x) = 2xF{ q \x) . (5) 

In the complete structure function, F^ (x) is weighted by the electric charge 
Q q of the quarks: 

Fi(x) = J2Q 2 q Fi q \x) . (6) 

In the following, we will only specify single-flavor structure functions and 
omit the superscript (q). For the mesons under consideration, the structure 
functions are identical for both flavors. The polarized structure function 
gi(x) gives the fraction of spin carried by quarks: 

9i{x) = ]^[q){x) - q\{x)] . (7) 

The structure function g 2 (x) does not have a parton model interpretation. 

The structure functions bi(x), b 2 (x), b 3 (x), and b±{x) are particular to 
spin-1 targets as the kinematical factors in eq. (2) involve the target po- 
larization vector to second order, a feature that does not occur for spin-| 
targets. In parton model language bi(x) and 62(2^) depend on the quark-spin 
averaged distributions q m = \(q™ + q™) only: 

6i(x) = q\x)-q\x) , (8) 
b 2 (x) = 2xb 1 (x) . (9) 

Thus bi(x) and b2(x) measure the difference in parton distributions of an 
m = 1 and m = target. This difference is due to the fact that in quantum 
field theory any Lorentz boost changes the particle content of a state. These 
changes differ for different spin orientations relative to the boost direction. 
For a model discussion of b x see e.g. [8]. 

The paper is organized as follows. In Sec. 2 we recall some results from 
the operator product expansion concerning the relevant operators and their 
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matrix elements. Sec. 3 describes the method we use to extract matrix ele- 
ments from three-point functions. The lattice implementation of this method 
is discussed in Sec. 4. Sec. 5 is devoted to questions of normalization and 
renormalization. In Sec. 6 we discuss our results, and Sec. 7 presents our con- 
clusions. Appendix A contains our conventions, Appendices B and C describe 
some technicalities. The reader who is not interested in the computational 
details may skip Sees. 3 and 4. 



2 Operators and moments of the structure 
functions 

The moments of structure functions can be related to the reduced matrix 
elements of certain local operators between pion or rho states. The local 
operators we consider are built from 7 matrices and covariant derivatives 
and have the general form in Minkowski space 

(M) "'"* = ^ T %-*/7"ifl"...iB'"fr , (10) 

of),.^ = _i_ G// ,^ 7 « 75i r i ...ir> / , , (id 

where ip is the quark field, and Gff> is a suitably chosen diagonal flavor 
matrix. The symmetrized derivative operators D are defined as 

D = D - D . (12) 

For a spin-0 particle, the momentum vector p is the only quantity the 
matrix element can depend on, and the reduced matrix element v n is defined 
by 

_ traces |p) = 2v n [p^ 1 ■ ■ -p^ - traces] . (13) 

The notation ■ ■ ■ /i n } denotes symmetrization in the indices fj,i, fi 2 , ■ ■ ■ , 
\i n . Expectation values of operators involving the 75 matrix vanish from 
symmetry considerations as the pion is a pseudoscalar particle. 

For a spin-1 particle, the structure of the matrix elements is more compli- 
cated due to the polarization degrees of freedom. Now both types of operators 
contribute 1 : 



(p,A|6W^ 1 -' 1 "> -traces|p,A) = 2S 



a n p^ ■■■pV" 



1 Note that we have corrected in (15) a misprint in [7]. 
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+ d n (e*^ 1 (p, A)e M2 (p, A) - ^p^) p^ ---p^ 



1 

3 j 



(p, \\0i M){ ^ n} - traces\p, \) 



2i 



m 1 - 



r n e^e;(p,X)e a (p,X)p T p^---p^ 



, (14) 



(15) 



S denotes symmetrization in the indices fj,i, . . . , ji n and removal of traces. 
The reduced matrix elements are a n for the polarization-averaged contribu- 
tion, d n for the polarized contributions, and r n for the operators involving 

75- 

By performing an operator product expansion of (1), reduced matrix ele- 
ments of local operators can be related to moments of the structure functions. 
We define the n-th moment of a function f(x) as 



.(/)= [ 1 x n - 1 f(x)dx . 
Jo 



(16) 



One then finds to leading order, which is twist two, the following represen- 
tation of the moments of the pion structure functions: 



2M n (Fi) 
M n _!(F 2 ) 



M 1 ) 



- v n 

for the rho structure functions one obtains [7]: 

2M n (Fx) = C^a n 



(17) 



M n _!(F 2 ) 
2M n (6x) 

M n _!(6 2 ) 
2M„(ffi) 



(18) 



where C^ fe ) = 1 + 0(a s ) are the Wilson coefficients of the operator product 
expansion. These relations hold for even n, except for the last one, which 
is valid for odd n. However, since we are calculating in the quenched ap- 
proximation, we are allowed to make use of these formulas for all n keeping 
in mind that our results can be meaningfully compared only with the non- 
singlet valence quark distribution. 

In the case of the pion, the moments of the quark distribution are given 

by 



(x 



n-l\ 



(19) 
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while for the rho they are related to the matrix elements a n : 

(x n " l ) = a n . (20) 

3 Three-point functions and matrix elements 

In order to calculate the reduced matrix elements on the lattice, we must 
calculate the expectation values of local operators of the form (10) and (11). 
To this end, we first need the connection between the Minkowski operators 
and those in Euclidean space. Defining Euclidean operators by 

OfX, n = Gtr ^rr M£? ) D^...D^ r , (21) 
C = 2"-i ^i? 75 D « • • • D > (22) 

we obtain the following relation to the operators in Minkowski space: 

Q(M)ni~nn _ ^_y*4+n 5 ^_^ni2 3 ; (23) 

where n 4 is the number of time-like indices, n i2 3 the number of spatial indices, 
and n 5 — 1 if the operator carries a 75 matrix. For our Euclidean conventions, 
see Appendix A. 

Lattice operators with the appropriate continuum behavior can be con- 
structed from the Wilson fermion fields by considering their symmetry prop- 
erties under the hypercubic group H(A) [9]. The operators we have chosen 
and their relation to the reduced matrix elements are listed in Appendix B. 

The required expectation values of our operators are extracted from ratios 
of two- and three-point functions. The three-point functions we consider are 
of the general form 

(rj(t^O(r)rj\0^) (24) 

where rj(t,p) is the sink operator for a particle moving with momentum p 
in time slice t, and rf(0,p) is the corresponding source at time slice t — 0. 
These operators are required to have the correct symmetry properties for the 
particles in question and their corresponding Hilbert space operators fj(p) 
should create the desired particles from the vacuum with nonzero amplitude. 
0(t) represents the operator O whose expectation value is to be calculated. 
For the pion we write 



(0\fKn;p)\n;p) = yJZ w , (25) 
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while for the rho there are three different particle states and correspondingly 
three different operators arranged in a vector r}i(p;t,p) that satisfy 



{0\fjj(p;p)\p;p, A) = Jz p €j{p,X) 



(26) 



up to lattice artifacts (see Appendix A for the definition of the polarization 
vectors e,). The correlation function for the rho depends on the polarization 
vectors: 

C jk = ( V ,(p;t,p)O(r)4(p;0,p)) . (27) 

In order to relate (27) to the matrix elements we are interested in, we 
express this correlation function in terms of traces involving the transfer 
matrix S: 



(r ]3 (p;t,p}O(r)4(p;0,p}} = 




S^ip^S^OS^p) 

s T ^ds^%{p-f)S%{p-^ 



T>t>r>0 
T > t > t > 



(28) 



Here T denotes the time extent of our lattice whose spacing is put equal to 
1. Then we insert a complete set of orthonormal eigenstates of the transfer 
matrix. If the time differences are chosen sufficiently large, we can restrict 
ourselves to the lowest contributing states \p, A), A labeling the three degen- 
erate polarization states of the rho. For the first case in (28) one obtains 



c jk = H( I%(p;p)Ip;p> a ) (p;p,Mo\p;p,\') {p;p,x\vk(p;p)\o)e 



-Et 



AA' 



= z pJ2 e j(P^)4(P^') {P',P,MO\p;p,X')e 



-Et 



AA' 

Z p m 2 T 3k e~ Et 



where Tj k is the matrix element between Cartesian states, 
T ik = \ J2 e j^ X ) e l(P^')(p-,p,\\0\p-,P,y) 



m 



AA' 



In the second case, one finds an additional sign factor: 



(2) 



f-1 



,n4+n 5 



Z p m T jk e 



-E(T—t) 



(29) 



(30) 



(31) 



To calculate all Cj k components for a given momentum would be ex- 
pensive in computer time. Choosing the momentum in 1-direction, we have 
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restricted ourselves to the components C 33 and C 32 . If (p;p, + \0\p;p, — ) = 
and (p;p,+\6\p;p,+) = (p;p, -\6\p;p, -), then 

(p;p,+\0\p;p,+)=T 33 , T 32 = , (32) 

whereas for (p;p, +\0\p;p, +) = -(p;p, -\0\p;p, -), 

(p;p,+\0\p;p,+) = -iT 32 , T 33 = . (33) 

The latter case is relevant to spin-dependent operators. To satisfy 

(p;p,+\6\p;p, -) = , 

it is sufficient that the operators commute with rotations in the plane trans- 
verse to p. This has also motivated our choice of operators. 

The factors that do not depend on the operator O can be eliminated by 
taking the ratio of (27) to another correlator, e.g. the two-point correlator 

j 

= E [(°I%(p;$Ip;p> a ) (p;p,Mvi(p;p)\Q)e~ Et 

A 

+(p-p : \\f) j (p-p)\0) (0|^(p;p)|p;p,A)e^ T -')] . (34) 
Using the relations (26) and (73), this reduces to 

C(t) = Z p (2m 2 + E 2 ) (e- Et + e^^) . (35) 

We therefore arrive at the following relation between the ratio of a three- 
to a two-point correlation function and the expectation value of the corre- 
sponding operator, valid for t 3> r 3> 0: 

(Vj(P^P)O(T)r,Up;0,p)) 



Rjk — 



Ei{vi(p;t,p)vi(p;Q,p)) 

1 e~ Et 



2 + E 2 /m 2 e- Et + e- E ( T -V Tjk ' (36) 

For T > r > t > 0, we get an analogous equation with the additional sign 
factor from (31) and with t replaced by T — t. For t = T/2, which is the 
choice in our numerical work, (36) gives 

Rjk = 2 + Ey m ^2 Tjk ■ (37) 
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The ratio may still depend on r due to contributions from the higher states 
neglected in (29). By searching for plateaus in the r-dependence, one can 
extract the value of the ratio with the smallest contamination from higher 
states. 

In the case of the pion, there is no polarization, and the relation (37) 
reduces to 

R= ^0\^ . (38) 



4 Evaluation of three-point functions on the 
lattice 

The actual form of the three-point correlator is given by 

( V F(t,p)O G (T)riF'(0,-d) . (39) 

Here we explicitly indicate the flavor matrices F, F', and G. ?7r(i,p) is a 
meson operator with momentum p at time t: 

V?(t,P)= E ^Fff^r^) (40) 

x:X4=t 

(a color, / flavor, a Dirac index) with a suitably chosen Dirac matrix T. In 
the case of the rho, T = 7^, while for the pion T = 75. A second meson 
operator is set at time slice with momentum —p, and F 1 = F' . The 
operator O g (t) has the general form 

G (r)= J2 G ff ,$%{z)J§{z,z>-x)^,p{z') (41) 

x,z,z':X4=t 

where J£p(z, z'\ x) is a matrix that represents the flavor, Dirac, and derivative 
structure of the corresponding local operator, x can be imagined as the 
"center of mass" of the operator while the sum over z and z 1 represents the 
derivative structure. 

Inserting these definitions, the correlation function is 

(^(t,p)O G (r) V ^(0,-p)) 

= ^3 E E E e ^^^v G ^ r ^r^ 

x:Xi=r y:y4=t z,z' 

x (J*(z, z'; x) r fa {y)r gP {y) ^(zWwM ^a'(0)^(0)> (42) 
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where V3 is the volume of a time slice. 

We integrate out the fermion fields in the quenched approximation and 
define 

(^*)^(y)> fc _ = 6 fr G%{U\x, y) (43) 

where G(U\x, y) is the fermion propagator in the gauge field configuration U 
(the U dependence will be indicated explicitly only when needed), and the av- 
erage is over fermion fields. There are six different contraction terms. In four 
of them, two operators at the same location are contracted. These fermion 
line disconnected contributions are proportional to trF, trF', and trG and 
vanish if these matrices are chosen traceless. However, this is in general 
impossible for trG, as we shall see below, and the omission of the corre- 
sponding contraction must be regarded an approximation, which is however 
consistent with quenching. We use this approximation for the same reason 
we use quenching: it is very hard to go beyond it. The remaining two terms 
are the fermion line connected contributions 

-v 3 E E 

x:X4=t y:y 4 =t z,z' 



(tr FGF') tr DC (TG(y, z)J(z, z' ■ x)G(z', 0)r'G(0, y)) 
+ (tr F'GF) tr DC (TG(y, 0)VG(0, z)J(z, z'; x)G(z', y)) ) (44) 



where the traces are over Dirac and color indices, and the average is over the 
gauge field alone. 

The two terms can be related to each other by means of the following 
relations 

G(x,yY = l5 G(y,x) l5 (45) 
7 5 r = sF^ 5 (46) 

75 r' = s'rS (47) 

7 5 J(z,/;a;) t 7 5 = sjJ(z',z;x) (48) 

where s, s', sj = ±1, and eq. (48) is valid only if the correspondig operator 
is suitably symmetrized in its space-time indices. For the pion, r = V = 75 
and thus s — s' — 1, while for the rho, T,V G {71,72,73} and s = s' = — 1. 
Then the correlation function reduces to 

-^3 E E e_W [( tr FGF')M(x, y) + ss'sj (tr F'GF) M(x, y)*] (49) 

x-.X4=T y:yi=t 
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with the basic single-flavor correlation function 

M(x,y) = J2( tr DcrG(y,z)J(z,z';x)G(z',0)r'G(0,y)) g . (50) 

z,z' 

Note that the calculation of this quantity on the lattice requires only two 
inversions of the fermion matrix, one at and one at y or z. 
Using the charge conjugation matrix, defined by 

rf = -C-\C , (51) 

and the relations 

G(U\x,y) = CGiU^xfC' 1 (52) 

CY T C- 1 = aT (53) 

CY' T C~ l = a'V (54) 

CJ(U*\z' ',z;x) T C- x = ajJ(U\z,z';x) (55) 

(where we explicitly denoted the dependence of J on the gauge field) with 
a, a', a j — ±1 one can further show that 

M(x,y)* = aa'ajss'sjM(x,y) . (56) 

We choose traceless matrices for F and F', 

'-(Si) • '-(!S) ■ (57) 

and therefore 

tr F'GF = Gn , tr FGF' = G 22 . (58) 

We finally arrive at the following expression relating the propagators 
M(x, y) to the three-point correlation function: 

-V 3 (G 11 + a J G 22 ) x ]T ]T e-**M(x,y) (59) 

X:X4=T y:y 4 =t 

For an operator with n derivatives, oj is (— l) ri+ ™ 5+1 ) where — 1 if the 
operator contains a 75 matrix, n5 = otherwise. Thus, for odd n + 715, G 
must not be traceless. 

The analogous expression for the two-point correlation function reads 

{r£(t,p)r£!(0,-p)) = 

= -V s E e^(tT DC G(x,0)T'G(0,x)T} g . (60) 
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5 Lattice and continuum operators 



Eqs. (36) and (38) relate the numerically computable ratios Rij and R to 
expectation values of Euclidean lattice operators. To connect them with the 
corresponding continuum Minkowski-space operators, we first introduce the 
continuum matrix element of the renormalized Euclidean operator O^ ont by 
the relation 

Z 6 (p\\0\p\) = ^ i- co ^A|6rVA> cont • (61) 

The factor 2E(p) is a consequence of the different normalization on the lattice 
and in the continuum: 

(p\\p'\') = Sptf,5 x ,\> 
c ^(p\\p'\r° nt = {2nf2E{p)5^5{p-p l ) , (62) 

and 2k comes from the definition of the Wilson fermion action on the lattice. 
Zq is the renormalization constant of the operator O. 

In the following, we shall use the renormalization constants calculated 
in one- loop lattice perturbation theory in the chiral limit [10]. They can be 
written in the form 

Zd = 1 " 16^ ° F ( ld lnM + Bd ~ B °o) ' (63) 
where Cp = 4/3, g denotes the bare coupling constant, and \i is the renor- 
malization scale. Note that here the lattice spacing a has been introduced 
explicitly. The finite contribution Bq is fixed in the momentum subtraction 
renormalization scheme, whereas B c ^ represents the contribution of the con- 
tinuum operator in the MS scheme with an anticommuting 75. Hence multi- 
plication by Zq leads from bare operators on the lattice to the corresponding 
renormalized (in the MS scheme) operators in the Euclidean continuum. For 
the renormalization scale fi we choose the inverse lattice spacing a -1 . Taking 
the physical rho mass of 770 MeV as input, we obtain from the lattice masses 
extrapolated to the chiral limit the value /i = 2.4 GeV. 



6 Results 

We have collected more than 500 independent configurations on a 32 x 16 3 
lattice at f3 — 6.0 with Wilson fermions and r — 1. Three different hop- 
ping parameters, k = 0.1515, 0.153, and 0.155 were used. They correspond 
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to quark masses of roughly 190, 130 and 70 MeV, respectively. As in [2], 
each gauge update consisted of a single 3-hit Metropolis sweep followed by 
16 overrelation sweeps. This cycle is repeated 50 times to generate a new 
configuration. The code was run on a Quadrics QH2 data-parallel computer. 
For completeness, the smearing technique - Jacobi smearing - we use to im- 
prove the overlap of the operator with the state is described in Appendix 
C. 

The calculational procedure is as follows: We calculate in each configura- 
tion the three-point functions (50) for a large set of operators as well as the 
pion and rho two-point functions. In Appendix B we list the operators we 
have actually studied. Those without 75 are labeled by the pion moments v n 
one can compute from them. The expectation value of such an operator Vn 
in the rho is a linear combination of a n and d n . The operators with 75 are 
labeled by the corresponding rho matrix elements r n . 

Using two values of the momentum, namely p = (0, 0, 0) and p = (||, 0, 0), 
we can check the continuum dispersion relation of the one-particle energies 
extracted from the two-point function. It is satisfied to better than 1%, and 
even for nonzero momentum we have a good projection on the ground state 
pion and rho. The particle masses we have used in our subsequent analysis 
are taken from Ref. [11]. They are collected in Table 1. 

For the computation of the three-point functions, the locations of the 
source and the sink are held fixed at and t — T/2 — 16. Placing the sink at 
T/2 allows us to search for a plateau equally well on both sides of the sink. 
In the case of the rho, we restrict ourselves to the 3-3 and 3-2 components. 

For the denominator of the ratios we employed two different procedures: 
first, we took the actual value of the propagator at midpoint, and second, we 
fitted the interior 24 points of the propagator to exponential functions and 
used the resulting midpoint value. The second case resulted in somewhat 
smaller errors at certain values of k and p l . We quote our results including 
the uncertainty from the former procedure. We also tried to use the conserved 
vector current, as proposed by [1], but this did not reduce our error margins. 

The ratios (36) and (38) are taken as a function of the operator insertion 
point r, and a fit to the central 7 points on each side that make up the 
plateau is made. The full covariance matrix is considered in estimating the 
error, thus taking correlations between neighboring points into account (in 
fact, only about 2 independent degrees of freedom out of 7 survived). Some 
example plots are shown in Figs. 1 and 2. 

In a few cases we have two operators for the same reduced continuum 
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Figure 1: The ratio R for the operators V 2a, Ov2b, V 2b at p ^ 0, O v s, and 
O v4: (left to right, top to bottom) for the pion at k — 0.153. The horizontal 
line is a fit to the central seven points on both sides. 
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Figure 2: The ratios -R33 and -R32 for the rho operators V 2 a , V 2b, O v ^, 
O r i a , and O r ib (left to right, top to bottom) at k — 0.153. 
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k = 0.1515 


k = 0.153 


k = 0.155 


k = k c = 0.15717(3) 


Pion 
Rho 


0.5033(4) 
0.5682(7) 


0.4221(4) 
0.5058(8) 


0.2966(5) 
0.4227(15) 




0.328(5) 



Table 1: Pion and rho masses in lattice units. 

matrix element, distinguished by the subscripts a and b. O v2a and O v2> b as 
well as O r 2,a and O r2 ,6 belong to different representations of the hypercubic 
group H(A). Hence the results extracted from operator a and operator b 
have to agree only in the continuum limit where the full 0(4) symmetry is 
hopefully restored, and a comparison of our results obtained for finite lattice 
spacing gives us some indication of the size of lattice artifacts. In the case of 
the operators associated to v 2> b, r i,b, and r 2ja we denote the results obtained 
with nonzero momentum by an additional subscript p. 

Applying (38), (61), (76), and (78) we have calculated estimates for the 
pion moments v n from the measured ratios R. The results are summarized 
in Fig. 3 and in Table 2. The agreement (within errors) of v 2)CL and v 2) b 
indicates that - at least in this case - lattice artifacts are not too large. 
Assuming a linear dependence on 1/k, i.e. on the bare quark mass, the values 
have been extrapolated to the chiral limit k — k c — 0.15717(3). Since the 
quark masses in our simulation are rather large (> 70MeV) we need this 
extrapolation in order to obtain numbers that can sensibly be compared 
with phenomenological valence quark distributions. Note, however, that the 
quark mass dependence of the results is not very pronounced. Only v 2> b shows 
a significant trend towards smaller values as the chiral limit is approached 
which is the expected behavior. 

We now come to the rho results (see Figs. 4, 5 and Table 2). By means 
of (32), (33), (61), and (76) we pass from the ratios (37) to matrix elements 
whose relation to a n , d n , and r n is listed in Eq. (79). The extrapolation to 
the chiral limit is performed as for the pion. In the case of the operators 
without 75, i.e. those labeled by v n , we encounter the problem that instead 
of one number (v n ) we have to extract the two quantities a n and d n from the 
matrix elements. 

Therefore we proceed as follows: The expectation value of O v2jb at p = 
gives us directly a 2 , and d 2 can then be calculated from O v2 , a and O v2t b at 
nonvanishing momentum. The expectation values of O v ^ and V 4, on the 
other hand, are proportional to d% and d 4 , respectively, if p = 0. With d% 
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<x> a <x> b <x> b . 



Figure 3: Estimates of the pion moments v n = (x^ 1 ) for a single flavor. For 
each matrix element, the results from the three different k values are shown 
versus 1/k with k decreasing (i.e. with the quark mass increasing) from left to 
right. The leftmost value is the chiral extrapolation obtained from a linear fit. 
The dotted lines give the free-field (heavy quark) limits. The dashed lines are 
phenomenological valence quark values from [3], evaluated at /i = 2.4 GeV. 



and d.4 computed from these matrix elements we use the corresponding results 
for p ^ to calculate 03 and a 4 . 

From the matrix elements of the operators with 75 we can easily extract 
r\, r 2 , and r 3 . The estimates r 2i(I , r 2jap , and r 2jfe for r 2 agree within the errors. 
Thus also in this case we do not observe significant discretization effects. 

7 Discussion 

We have calculated the lowest three moments of the structure functions of the 
pion and the rho meson, restricting ourselves to the leading twist-2 operators 
in the operator product expansion. 
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<x> 



2 

<X > 



3 

<X > 




d2 



■i — $- 



d3 




d4 



Figure 4: Estimates for the rho moments a n = (x n x ) and d n . The presenta- 
tion of the data is the same as in Fig. 3. 
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k = 0.1515 


k = 0.153 


k = 0.155 


k = k c = 0.15717 


PION 


V2,a = (x) a 


0.301(20) 


0.294(28) 


0.290(71) 


0.279(83) 


V 2 ,b = {x) b 


0.3239(70) 


0.3150(71) 


0.2910(75) 


0.273(12) 


V2,bp = {x)bp 


0.319(23) 


0.316(33) 


0.325(84) 


0.318(98) 


v 3 = (x 2 ) 


0.1222(83) 


0.116(12) 


0.117(31) 


0.107(35) 


v A = (x 3 ) 


0.0619(45) 


0.0580(65) 


0.054(18) 


0.048(20) 


RHO 


a 2 = (x) 


0.3555(80) 


0.3531(93) 


0.340(14) 


0.334(21) 


a 3 = (x 2 ) 


0.1398(93) 


0.144(14) 


0.182(48) 


0.174(47) 


a 4 = (x 3 ) 


0.0725(72) 


0.069(12) 


0.074(41) 


0.066(39) 


d 2 


0.107(52) 


0.128(75) 


0.29(20) 


0.29(23) 


d 3 


0.0145(32) 


0.0135(49) 


-0.002(14) 


0.001(15) 




0.0109(100) 


0.004(17) 


0.007(62) 


-0.009(58) 


n, a 


0.709(56) 


0.715(97) 


0.42(34) 


0.57(32) 


r hb 


0.721(17) 


0.702(20) 


0.627(32) 


0.590(46) 


r hbp 


0.680(56) 


0.62(13) 


0.32(44) 


0.33(42) 


T 2 ,a 


0.2743(62) 


0.2631(70) 


0.231(12) 


0.212(17) 




0.257(17) 


0.243(25) 


0.216(69) 


0.198(76) 


r 2 ,b 


0.242(20) 


0.232(30) 


0.210(89) 


0.199(95) 


r 3 


0.1067(71) 


0.099(11) 


0.087(33) 


0.077(34) 



Table 2: Result overview for a single flavor. The numbers refer to the MS 
scheme with a renormalization scale /i m 2.4 GeV. The last column gives the 
result of the extrapolation to the chiral limit. 
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Figure 5: Estimates for the rho moments r n . The presentation of the data is 
the same as in Fig. 3. 



For the pion, we can compare our numbers in the chiral limit with the 
experimental data [3]. Our result for (x) is larger than phenomenology sug- 
gests. This is to be expected as our quenched lattice calculation does not 
contain any sea quarks and the valence quarks will therefore carry more of the 
momentum. The results for (x 2 ) and (x 3 ), on the other hand, are consistent 
with the phenomenological numbers. Our results also agree with the early 
lattice calculations of Martinelli and Sachrajda [1] as well as with various 
model calculations [12]. 

The unpolarized rho structure function looks very similar to the pion 
structure function, at least for the quark masses that we have considered. In 
the pion the quarks carry about 60% of the total momentum, while in the 
rho they carry about 70% at the smallest quark mass. The higher moments 
are in agreement with each other within the error bars. Thus the assump- 
tion F[(x) ~ F[(x) often used in phenomenological estimates may well be 
justified. 

The lowest moment r\ of the polarized structure function g\ indicates 
that the valence quarks carry about 60% of the total spin of the rho. For 
comparison, a similar quenched calculation for the nucleon gave a quark spin 
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fraction of about the same value [2], which is reduced to 18% by sea quark 
contributions [13]. It is very likely that the same will also happen here. 

The structure functions 61,62 measure the difference in quark distribu- 
tions of a (spin projected) m — 1 and m = rho meson. If the quarks 
were in a relative s-wave state in the infinite momentum frame, we would 
expect 61, b 2 to be zero. The lowest moment d 2 turns out to be positive and 
surprisingly large on the scale of a 2 , albeit with large statistical errors. Per- 
haps this indicates that the valence quarks have a substantial orbital angular 
momentum. This could also explain a relatively small quark spin fraction. 
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A Conventions 

The Minkowski space metric has the signature (1, — 1, — 1, — 1). Minkowski 
and Euclidean components are related by 

X4 = ix (M)0 = ix (M) ^ Xj=x (M) J = _ X (M) ^ (64) 

where j refers to spatial indices. Unless explicitly mentioned, we label 
Minkowski-space variables by an (M). 

Covariant derivatives are defined in Minkowski space as 

D {M)n = g(% _ i g A^ (65) 
and are related to their Euclidean counterparts as follows: 

iD (M)0 = ^ iD (M)j = _ iD , ^ (66) 

similarly for the gamma matrices: 

7 (M) ° = 74 , 7 {M)j = i7, • (67) 
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The 75 matrix is defined as 

7 ^ M) = i 7 ( M )° 7 ( M ) 1 7 ( M ) 2 7 (M)3 

75 = 7i727374 = ~75 M) • (68) 

The momentum of the particles is chosen in the 1-direction, p = (p, 0, 0). 
Polarization vectors for vector particles satisfy 

P^fip, A) = , e^(p, A)eW(p, A') = -m 2 5 XX , (69) 

(A = ±, 0) and have the explicit form 

e {M Hv, A) = ( P -^, e A + -^^p) (70) 
\ m m\m + E) J 

with the basis vectors 

e± = T-^(0,l,±i) (71) 
e = m (1,0,0) . (72) 

They satisfy the completeness relation 

E ^V, A)6^(p, A) = m 2 (<5 y + -L m ) . (73) 

Note that in Euclidean space A) = e^ M ^(p, A). 

B Operators 

On the lattice, the choice of the operators to look at is a nontrivial matter, 
because the discretization reduces the symmetry group of (Euclidean) space- 
time from 0(4) to the hypercubic group H(A) C 0(4). Hence the lattice 
operators have to be classified according to H(4) and one should choose op- 
erators belonging to a definite irreducible representation of H(4). Since H(A) 
is a finite group, the restrictions imposed by symmetry are less stringent than 
in the continuum and the possibilities for mixing increase. Whereas mixing 
with operators of the same dimension is supposed to be treatable by pertur- 
bation theory, the mixing coefficients for lower-dimensional operators have to 
be calculated nonperturbatively. Hence one would like to avoid mixing with 
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Operator 


z 6 


"X 

O v 2,a 


0.989196 


O v 2,b 


0.978369 


0,3 


1.102397 


V 4 


1.229911 



Operator 


Z 6 


^ 

O r l,a, O r ifi 
O r 2,a 
O r 2,b 

o r3 


0.866625 
0.997086 
0.998587 
1.108573 



Table 3: Renormalization constants Z. 

lower-dimensional operators whenever possible. On the other hand, as the 
spin grows, operators with no mixing at all require more and more nonva- 
nishing momentum components in the calculation of their forward hadronic 
matrix elements, which makes their Monte Carlo evaluation increasingly dif- 
ficult. So some kind of compromise is needed. 

Due to our use of the quenched approximation purely gluonic operators 
cannot mix with two-quark operators and we may restrict ourselves to the 
latter. Guided by their H(A) classification given in Ref. [9] we have chosen 
the following operators in Euclidean space: 



O v 2,a 
O v 2,b 

0,3 
0„4 
O r l,a 
O r l,b 
O r 2,a 

O r 2,b 

o r3 



o 



{41} 



6 44 - \ (On + O22 + O33) , 



0{H4} 



o 



{1122} 



1 

~ 2 

+ 6 



o 



{224} 



o 



{334} 



{3344} 



+ 



{1133} 



+ 



{2244} 



-20 



{1144} 



-20 



{2233} 



of 



oi 

6 5 



{41} 



4 5 4 - l (Oil + Ol2 + 3 5 3 ) , 



^{114} 



7 {22 1} + Oi'; !:i |}) 



(74) 



For V2 and r2 we have two operators, which belong to the same 0(4) multiplet 
in the continuum limit but transform according to inequivalent representa- 
tions of H(4). Hence their matrix elements provide a test for the restoration 
of 0(4) symmetry. The renormalization constants for these operators in the 
MS scheme are listed in Table 3. 
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Concerning the mixing properties a few remarks are in order. Mixing 
with operators of equal or lower dimension is excluded for the operators 
O v2 , a , O v2 ,b, O rlta , O r i tb , O r2 , a , O r2 ,b- The case of the operator 0^3, for 
which there are two further operators with the same dimension and the same 
transformation behavior, is discussed in ref. [9]. The operators 0„4 and O r 3, 
on the other hand, could in principle mix not only with operators of the same 
dimension but also with an operator of one dimension less and different chiral 
properties. It is of the type 

^<V75 D m D H2--- D ^ , (75) 

where n = 2 in the case of V 4, and n — 1 for O r3 . 

Our analysis ignores mixing completely. This seems to be well justified for 
V 3. Here a perturbative calculation gives a rather small mixing coefficient 
for one of the mixing operators, whereas the other candidate for mixing does 
not appear at all in a 1-loop calculation, because its Born term vanishes in 
forward matrix elements. The same is true for all operators of dimension less 
than or equal to six which transform identically to O v ^: Their Born term 
vanishes in forward matrix elements, hence they do not show up in a 1-loop 
calculation. In the case of O r3) however, the mixing is already visible at the 
1-loop level. The results for t> 4 and r 3 have therefore to be considered with 
some caution. 

The corresponding Minkowski operators are found by applying eq. (23). 
Defining the Minkowski analogs of our Euclidean operators by 

A - \A(.M) A _ A(M) 

A _ A(M) A _ A(M) 

V V 3 — —V v3 , U V 4 — U v4 , 

Orl,a = -Offl , Orl,b = iO™ » ^ 

A _ ;A(M) A A(M) 

<^r2,a — — "-A-2,a j ^r2,b ~ U r 2,b > 

n n (M) 

Vr3 — U r3 , 

we have 

AW = A(M )m > 



= ()(M)00 + 1 ^(M)ll + q(M)22 + q(M)33^ 

o { ™ ] = d {M){lw} - - (p {M){220} + d (M){330} ) 
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_ _ q(M){1122} _|_ £j(M){3300} _ q(M){1133} + £)(M){2200} 

- 2d (JVf){1100} + 20 (M){2233} 

A(M) _ A(M)0 

<-'rl,a — ^5 J 

A(M) A(M)1 

°Vl,6 — U 5 ) 

A(M) _ A(M){01} 

oSS = of /)0 V i (of )n + of )22 + of )33 ) , 

( f ) = 6f ){110} - - (6f ){220} + 6f ){330} ) (77) 

We can now use (13), (14), and (15) to calculate the expectation values of 
the Minkowski space operators. For the pion, one obtains with p 2 = p 3 = 

<oS2> = WV , 

<oS3> = 2v 2 ((py + 1 -( P r) , 
(of)) = 2^)y , 

(Of } ) = -^(P 1 ) 2 ^ ) 2 • (78) 
For the rho with polarization A = ±, one finds 

(Of!) = 2(a 2 -\d 2 )pY , 

(OS) = 2(a 2 -^ 2 )((p°) 2 + i(^) 2 ) + ^ 2 , 

(O22) p i=o = 2a 2 m 2 , 

(Of } ) = 2(a 3 -^)b 1 )V-^V , 

(Of)) = -4 (a 4 - ld 4 ) (^V) 2 + Vrf 4 , 

(Of j) = ±2^ , 
(Of)) = ±2r lP ° , 

(Ofj) = ±r 2 ((p°) 2 + (pr) , 
(Of)) = ±j^V , 

(Of)) = ±^ 1 ((p 1 ) 2 + 2(p ) 2 ) . (79) 
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C Smearing 



The method we use for smearing is to smear the quark in a plane £4 = t [14]: 

s r fa (x, t) = e s # a6 ( f > ^ w)M t) , (so) 

where the kernel H is chosen to have the correct gauge transformation prop- 
erties and is diagonal in spin space. S is the smearing label. So for example 
for no smearing S = L = local, then L H ab (x,y;U,t) = 5 ab 5^. S H is also 
taken as Hermitian: 

s H ba (y, x- U, t)* = s H ab (x, y- U, t) . (81) 

Also the smeared anti-quark is defined as 

s r fa (x , t) = £ ty a (y, t)H ba (y, x- U, t) . (82) 
y 

Note that we can choose different smearing for the quark and anti-quark. 
Thus for a smeared meson operator we have 

S ' S v(t,P) = ^F ffl e-^ s 'i,%{x^)T aP s r r ^t) (83) 

with the appropriate correlation function 

C s ' s (t,p;t ) = ( S ' S 'v(t,P) SS v(to,-p)) (84) 

so that S and S' is smearing at the source, sink respectively. 
The smeared quark propagator is defined by 

s ' s G f a f(x, y; U)5 ff , = ( S > / a a (x)^^(l/))fe rm io„s (85) 

(LLq = Qy g Q j n meson correlation functions we can simply replace G with 
S ' S G to allow for smearing. 

The smeared quark propagators are found sequentially: 

• Generate the smeared source S S from a point source at (x ,t ) and so 
with S$ a (xt) = 5^ Q 5 ttQ 5 aa "5 aao , 

s SlZ° (x, t; x , t ) = s H aa ° (x, x ; U, t)S tto 5 aao . (86) 
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• Find LS G, by solving M LS G = S S. We thus have 

LS G f a f(x, t; x , t ) = ]T G f a f(x , t; y, t ) s H b ' b (y, x ; U, t ) . (87) 

y 

• From LS G we generate S ' S G by applying S 'H: 

s ' s G f a f(x, t- x , t ) = £ s 'H aa \x , y; C/, t) L5 Gg b (y, t; f , <b) ■ (88) 

Note that this step can be expensive (in CPU time) in comparison to 
eq. (86) as we must smear on every x 4 = t plane. 

Practically we shall use Jacobi smearing (as advocated mainly by [15]). 
This is given by 

K{x, t; xt) s S(x', t) = S (x, t) , (89) 

x' 

where So is the original point source. Here 

K — 1 — k s D s (90) 
and D s is a covariant derivative in the x± — t plane, viz. 

D%p(x, t; y, t) = 5 aP 2 [Uf(x, t)5 3+ ^ + Uj ab (x - 7, t)fe_,J . (91) 
i=i 

Hence we need H = K~ x . Rather than performing this inversion completely 
we Jacobi iterate N s times, so 

s S (n \x,t) = S {x,t) + K s D s s S (n ' 1 \x,t) n = l,2,... (92) 

with s S(°\x,t) = S {x,t). 

We thus have two parameters, k s , N s at our disposal. k s controls the 
coarseness of the iteration, while increasing N s increases the size of the 
smeared object roughly like a random walk. Physically we wish to smear 
until our source is about the size of the meson. A suitable measure of the 
(rms) radius is given by 

2 



r 



2 



Y,x{x-xof s S(x,t ;x ,t ) 



(93) 

\ s S(x,t -x ,t )\ 

(Note that on a periodic lattice, (x — Xq) 2 is taken as the minimum distance 
from x to x .) Explicitly for (3 = 6.0 we have chosen k s = 0.21, iV s = 50. 
This gives ra of about 3.5a > 0.5fm which corresponds roughly to the hadron 
radius. 
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